PROCEEDINGS 

OP SCIENCE 



Charmed meson spectroscopy on the lattice 



Daniel Mohler*", R. M. Woloshyn'' 

" TRIUMF,4004 Wesbrook Mall, Vancouver, BC V6T2A3, Canada 
E-mail: jnohlergtriumf . ca rwww@triumf . ca 



We present a calculation of the spectrum of low-lying D and mesons calculated on 2+1 fla- 
vor configurations with Clover- Wilson quarks generated by the PACS-CS collaboration. S- and 
P-wave states are explored for pion masses down to 156MeV, including some excited states. As a 
benchmark, a calculation of the low-lying charmonium spectrum is performed. While the charmo- 
nium results agree favorably with experiment, noticeable differences remain for charmed mesons 
in channels where resonances close to multi-particle thresholds exist. 
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1. Introduction 

While the charmonium spectrum below the DD and D*D thresholds is well described in the 
quark model [jl]], the spectrum of charmed mesons contains states for which model expectations 
did not hold. Of special interest are the charmed strange mesons Z)*q(2317) and Ds\{2A60) which 
would form a mass-degenerate pair with = ^ (with P parity and j the total angular momentum 
of the light quark) in the heavy quark limit. While model expectations placed these states above 
the DK and D*K thresholds respectively, the experimental states have masses below threshold and, 
as a consequence, are very narrow. This experimental fact sparked speculations about the nature 
of these states. Lattice QCD is ideally suited to determine the low-lying spectrum from first prin- 
ciples. Early lattice simulations employed the quenched approximation or very heavy sea quarks 
and most simulations observed a large discrepancy with regard to experiment. In the following, 
results obtained on dynamical gauge configurations with pion masses as light as 156MeV in a large 
box of size 2.9fm are presented. Details have previously been published in [^. In section ^ we 
describe the details of our simulation and in section || our main results for the low lying spectra of 
charmonium and charmed mesons are presented. 

2. Calculational setup 

The Clover- Wilson gauge configurations made available by the PACS-CS collaboration 
were used. In the simulation the strange quark is fixed close to its physical mass and pions made 
from light (up and down) quarks range from 702MeV to I56MeV across the different ensembles. 
The lattice spacing a = 0.0907(13) has been determined in and the lattices are of extent 32-' x 
64. Table [l| lists the run parameters and the number of configurations used in this study. 



Ensemble 


Jh) 






#configs D/D, 


1 


1.52617 


0.13700 


0.13640 


200/200 


2 


1.52493 


0.13727 


0.13640 


-/200 


3 


1.52381 


0.13754 


0.13640 


200/200 


4 


1.52327 


0.13754 


0.13660 


-/200 


5 


1.52326 


0.13770 


0.13640 


200/348 


6 


1.52264 


0.13781 


0.13640 


198/198 



Table 1: Run parameters for the PACS-CS lattices All gauge configurations have been generated with 
the inverse gauge coupling j3 = 1.90 and the light quark clover coefficient ci2 = 1.715. The quantity ciw 
denotes the heavy quark clover coefficient used for the charm valence quarks. 

For the heavy charm quark the Fermilab method [Q, ^ is employed in a fashion very similar 
to [^. The heavy quark hopping parameter Kc has been tuned for the spin-averaged kinetic mass of 
the charmed-strange IS states to assume its physical value. The heavy quark clover coefficient is 
set to its tadpole improved value 4, , where uq is the average link as determined from the plaquette. 
For a detailed description of our tuning please refer to 
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Figure 1: Left panel: IS hyperfine splitting compared to the physical splitting. Right panel: Splitting 
between the spin-averaged S- and P-wave states. Lattice errors are statistical only. 



The low-lying spectrum is extracted using the variational method For each channel a 

correlation matrix C{t) with a number of lattice interpolating fields of the desired quantum numbers 
(or J^^ for charmonium) is constructed 

Cit)ij = £e-'^" (O|0,-|n) {n\o]\0) . (2.1) 



On each time slice the generalized eigenvalue problem (GEVP) is solved 

C(Ov)^W =A«(f)C(?o)v>(*\ (2-2) 
A«(0-e-'^*(l + ^(e-'^*)). 

Ordering the eigenvalues by magnitude or by their eigenvectors, the ground state can be extracted 
from the asymptotic behavior of the largest eigenvalue, the first excited state from the second largest 
eigenvalue and so on. Details of our interpolator basis, which contains both Jacobi-smeared [Io| ] 



Gaussian-shape sources and derivative sources [11, 12] can be found in p 



The quark propagators are calculated for sources at random spatial points on eight different 
source time slices. For the calculation of the light-quark propagators the dfl_sap_gcr inverter from 
Liischer's DD-HMC package 14] is used. For the charm quark propagators the corresponding 



inverter without deflation is employed. For the large number of sources and light dynamical quarks 
the use of a deflation inverter was a crucial ingredient. 

3. Results 

The low-lying charmonium spectrum below multi-meson thresholds is an excellent test for our 
setup, in particular for the tuning of the charm quark hopping parameter which has been determined 
using information from the charmed-strange spectrum. Furthermore, all charmonium states below 
the DD* and DD threshold are well determined and commonly believed to be regular qq states. 
To illustrate the quality of our results Figure [| shows selected results. The left panel shows the 
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Figure 2: Mass splittings in the charmonium spectrum compared to the spin-averaged ground state mass 
{M-q^ + 3M/>j/)/4. The errors have been obtained by combining the statistical and scale setting uncertainties 
in quadrature. 

chiral behavior of the charmonium hyperfine splitting. While the statistical errors are tiny, strong 
discretization effects are expected for spin-dependent quantities. It is therefore no surprise that 
we underestimate the hyperfine splitting by almost 20MeV at this lattice spacing. The right panel 
shows the splitting between the spin-averaged S- and P-wave states. While the spin dependent 
splittings are reproduced rather poorly the spin-independent splitting is completely consistent with 
the experimental value. 



Mass difference 


Our results [MeV] 


Experiment [MeV] 


IS hyperfine 


97.8 ±0.5 ±1.4 


116.6±1.2 


IP spin-orbit 


37.5 ±2.4 ±0.5 


46.6 ±0.1 


IP tensor 


10.44±1.13±0.15 


16.25 ±0.07 


2S hyperfine 


48±18±1 


49±4 



Table 2: Spin dependent mass splitting in the Charmonium spectrum. 



An overview of the results for low-lying charmonium states can be found in Figure 0, where 
results from the lowest pion mass are plotted. Overall there is good qualitative agreement with 
the experimental spectrum. In addition Table ^ lists the spin-dependent splittings, which are very 
sensitive to discretization effects. 

Turning to the spectrum of charmed mesons. Figure ^ illustrates the importance of interpolator 
mixing in the 1+ channel, where two low-lying states are expected. In this case the variational basis 
with a combination of interpolators corresponding to positive and negative charge conjugation in 
the mass-degenerate case was crucial. This mixing enhances the mass-splitting between the two 
states, although our final results (see Figure ^ for this splitting are still to small. The left panel of 
Figure^ shows the spin-averaged splitting between the IS and 2S states in the spectrum. So far 
only one of these has been unambiguously determined in experiment. The right panel shows the 
hyperfine splitting between the 2S states. To our knowledge this is the first lattice determination of 
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Figure 3: Left panel: Selected single diagonal correlators for Dj in the 1+ channel. Right panel: Lowest 
two energy levels from a 4 x 4 matrix analysis. The data is from ensemble 6 
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Figure 4: Left panel: Mass splitting between the spin averaged 2S and IS states in the Ds spectrum. Right 
panel: Hyperfine splitting for the 2S states. Lattice errors are statistical only. 



the 2S states in the spectrum. 

Our spectrum results for the charmed mesons are summarized in Figure |5[ The data points are 
from the simulation with lightest pion mass. Overall we obtain reasonable results for the hyperfine 
splittings and for the pairs of states corresponding to the multiplet with = | in the heavy quark 
limit. The dynamical calculation with light sea quarks improves the overall agreement of the spec- 
trum with experiment. Nevertheless substantial differences remain for the doublets corresponding 
to 7'^ = J in the heavy quark limit, especially in the D., spectrum. For these states the nearby DK 
and D^'K thresholds may play an important role. Figure |6| compares the experimental masses and 
scattering thresholds with the corresponding results obtained in the lattice simulation. The scat- 
tering thresholds in our simulation are slightly unphysical^ . At the smallest pion mass the energy 

'The main discrepancy arises due to the shghtly unphysical strange quark mass. In addition the light quark masses 
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Figure 5: Mass splittings in the D and Dg meson spectrum compared to the respective spin-averaged ground 
state mass. The errors have been obtained by combining the statistical and scale setting uncertainties in 
quadrature. 



level determined in our simulation coincides with the scattering threshold. While this suggests that 
we may see the scattering level as the ground state in our simulation, the eigenvectors extracted 
from the GEVP suggest that we observe the same state for all ensembles used. We hope to revisit 
this issue in future simulations by including the relevant scattering levels in our interpolator basis. 
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Figure 6: Measured energy levels for the Z)*q (left panel) and D^x (right panel) ground states (black circles) 
compared to experimental states (magenta circles). All masses are plotted with respect to the spin-averaged 
Ds ground state. The plus signs denote the DK and D*K scattering levels on the lattice (black) and in nature 
(red). At our lowest pion mass the artificially heavy scattering states are very close to the measured ground 
state energy. 



are also slightly larger than in nature. 
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